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A Brownian particle having a fluctuating mass 
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We focus on the dynamics of a Brownian particle whose mass fluctuates. First we show that 
the behaviour is similar to that of a Brownian particle moving in a fluctuating medium, as studied 
by Beck [Phys. Rev. Lett. 87 (2001) 180601]. By performing numerical simulations of the 
Langevin equation, we check the theoretical predictions derived in the adiabatic limit, and study 
deviations outside this limit. We compare the mass velocity distribution with truncated Tsallis 
distributions [J. Stat. Phys. 52 (1988) 479] and find excellent agreement if the masses are chi- 
squared distributed. We also consider the diffusion of the Brownian particle by studying a Bernoulli 
random walk with fluctuating walk length in one dimension. We observe the time dependence of 
the position distribution kurtosis and find interesting behaviours. We point out a few physical cases 
where the mass fluctuation problem could be encountered as a first approximation for agglomeration- 
fracture non equilibrium processes. 

PACS numbers: 05.40.Jc, 05.20.-y, 89.75.-k 



I. INTRODUCTION 

It is well known that, while the Maxwell-Boltzmann 
distribution takes place in any system at equilibrium, 
non-equilibrium systems present in general qualitative 
and quantitative deviations from the former. A case of 
particular interest is that of distributions characterized 
by a power law tail, and therefore by an over-population 
for high energy as compared to the Gaussian. Indeed, 
such behaviours occur in a large number of physical sys- 
tems, going from self-organized media to granular gases, 
and may have striking consequences due to the large 
probability of extreme/rare events. 

The ubiquity of these fat tail distributions in nature 
has motivated several attempts in the literature to con- 
struct a general formalism for their description, one of the 
most recent being the Tsallis thermodynamic formalism 
P, 0] . The latter is based on a proper extremization of 
the non extensive entropy 
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that leads to generalized canonical distributions, often 
called Tsallis distributions 
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In this expression, x denotes the state of the system, and 
e x q is the ^-exponential function defined by 
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This definition implies that ef = e x , and that the en- 
ergy is canonically distributed in the classical limit q — ► 1 , 
as expected. Let us also note that this formalism draws 
a direct parallelism with the equilibrium theory, where 
(3 plays the role of the inverse of a temperature, and Z 
that of a partition function. 

Tsallis distributions have been observed in numerous 
fields [3 but their fundamental origin is still debated due 
to the lack of simple model and exact treatment justifying 
the formalism. It is therefore important to study simple 
statistical models in order to show in which context Tsal- 
lis statistics apply. There are several microscopic ways 
to justify Tsallis statistics. One of them has been in- 
troduced some years ago by Beck 1] . Another has been 
recently introduced by ThurnerQ. Beck theory was first 
initiated by considering the case of a Brownian particle 
moving in a specially thermally fluctuating medium, i.e. 
the inverse temperature is a chi-squared distributed 
random variable. Yet, Thurner has shown that one 
can derive Tsallis distributions in a general way without 
Beck " chi-square assumption" . 

In view of the above considerations, the motion of 
a Brownian particle whose mass is a random variable 
seems to be a paradigmatic example. Moreover many 
physical cases are concerned by such situations. A sharp 
mass or volume variation of entities can be encountered 

^3 

in many non-equilibrium cases: ion-ion reaction[5l |fj, |]|, 
electrodepositionQ, granular flow[j| fToL [ill formation 
of planets through dust aggregation^^ Il3l IbH ] , film 
depositi on!!!!!, traffic jams 

|lfl LL2P ]> and even stock 
markets |18|, |l9j (in which the volume of exchanged shares 
fluctuates and the price undergoes some random walk). 
It should be noticed that a priori fluctuating mass prob- 
lems differ from temperature fluctuations. Two masses 
can be addded to each other; this is hardly the case for 
temperatures. An expose of such a generalized Brow- 
nian motion and the distinction between masses and 
temperatures will be emphasized in Sect. II and Sect. III. 
By supposing two time scales, i.e. assuming that relax- 



2 



ation processes for the particle are faster than the char- 
acteristic times for the mass fluctuation, it will be shown 
that its asymptotic velocity distribution is in general non- 
Maxwellian. Moreover, for some choice of the mass prob- 
ability distribution, the Brownian particle velocities are 
Tsallis distributed. We verify this result by performing 
simulations of the corresponding Langevin equation. 

Moreover, we also consider the case when the mass fluc- 
tuation is not adiabadically slow. We show in Sect. IV 
that the mass dependence of the relaxation rate may have 
non-negligible consequences on the velocity distribution 
of the Brownian particle. In order to fit the resulting 
distributions, and to describe the deviations from Tsallis 
statistics, we introduce truncated Tsallis distributions. 
Finally, in Sect. V, we study the diffusive properties of 
the Brownian particle. To do so, we model the motion 
of the particle by a random walk with time-dependent 
jump probabilities, associated with the mass fluctuations 
of the particle. The model exhibits standard diffusion, 
i.e. < x 2 >~ t, but the shape of the scaling position dis- 
tribution is anomalous and may exhibit quasi-stationary 
features. Sect. VI serves as a summary and conclusions. 



g(f3) — ► 6(/3 — Pe)- In contrast, different choices for the 
statistics of f3 may lead to a large variety of probability 
distributions for the Brownian particle velocity. 

Several forms for g{/3) have been studied in the liter- 
ature _2 1] , but one functional family of g(/3) is particu- 
larly interesting. Indeed, the generalized Langevin model 
Eq. (0} generates Tsallis statistics for the velocities of the 
Brownian particle if (5 is a chi-squared random variable 

°W = W® ^ ^ (6) 

where b and c are positive real parameters which account 
for the average and the variance of j3. Let us stress that 
a chi-squared distribution derives from the summation 
of squared Gaussian random variables X i: 8 — J2i=i ^-li 
where the Xi are independent, and < Xi >— 0. By intro- 
ducing Eq.@ into Eq.JSJ), it is straightforward to show 
that the velocity distribution of the particle is Eq.©, if 
one identifies c = —K- and be = 8 . 

For completeness, let us also mention a study of Eq.(4) 
when A fluctuates 221. 



II. GENERALIZING THE BROWNIAN 
MOTION TO ONE IN A FLUCTUATING 
MEDIUM 

For setting the framework of the present study, let us 
recall one microscopic way to justify Tsallis statistics [lj. 
Beck theory, later called Superstatistics, was first initi- 
ated by considering the case of a Brownian particle mov- 
ing in a fluctuating medium, i.e. an ensemble of macro- 
scopic particles evolving according to Langevin dynamics 

md t v = -Xv + aL(t) (4) 

where A is the friction coefficient, a describes the strength 
of the noise, and L(t) is a Gaussian white noise Con- 
trary to the classical Brownian motion, however, one may 
consider that the features of the medium may fluctuate 
temporally and/or spatially, namely the quantity /3 = 
i.e. a quantity which plays the role of the inverse of tem- 
perature, changes temporally on a time scale r, or on the 
spatial scale L; see also [2jj. E.g. in his original paper, 
Beck assumed that this quantity fluctuates adiabatically 
slowly, namely that the time scale r is much larger than 
the relaxation time for reaching local equilibrium. In 
that case, the stationary solution of the non-equilibrium 
system consists in Boltzmann factors e~@ x that are av- 
eraged over the various fluctuating inverse temperatures 

fBeck(x) = ^fd(3 g((3) e-P* (5) 

where if is a normalizing constant, and g([3) is the 
probability distribution of 8. Let us stress that or- 
dinary statistical mechanics are recovered in the limit 



III. FLUCTUATING MASS 

Let us now consider the diffusive properties of a macro- 
scopic cluster such as one arising in granular media 0, IT(1 ] 
or in traffic 0|. Such media are composed by a large 
number of macroscopic particles. Due to their inelastic 
interactions, the systems are composed of very dense re- 
gions evolving along more dilute ones. In general, there 
is a continuous exchange of particles between the dense 
cluster and the dilute region, so that the total mass of 
the macroscopic entity is not conserved. As a first ap- 
proximation, we have thus considered the simplest ap- 
proximation for this dynamics, namely the cluster is 
one Brownian-like particle whose mass fluctuates in the 
course of time. To mimic this effect, we have assumed 
that (i) the distribution of masses is a priori given by 
g{m), and (ii) the mass of the cluster fluctuates with a 
characteristic time r. By definition, this model evolves 
according to the Langevin equation Eq. Q , where m is 
now the random variable. 

Given some realization of the random mass, say m — 
rriR, one easily checks that the velocity distribution of 
the particle converges toward the distribution 



f B {v) -+ -j=e , (7) 

This relaxation process takes place over a time scale 
in ~ Therefore, if the separation of time scales 

tii << r applies, the asymptotic velocity distribution of 
the cluster is given by 

fsiy) = J dmg{m)^e-^ . (8) 
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FIG. 1: Time evolution of the < v 2 > for three values of r, 
characterised by s = ~& = [0.1, 1, 10]. The possible masses in 
the system are 1 and 10. Consequently, the < v 2 > for each 
species are respectively 1.0 and 0.1. Moreover, in the slow 
limit (r >> tu), the asymptotic velocity fluctuation is, on 
average, 0.55. When r is small, most of the particle velocities 
are distributed like the light ones, and their energy is closer 
to 1.0 than to 0.55. 



Consequently, this leads to a Tsallis distribution if the 
masses are chi-squared distributed. In that sense there is 
a direct correspondence between Beck approach and ours. 
However there is more to see in the latter case because it 
justifies Tsallis non extensive entropy approach in a more 
mechanistic way. Subsequently two basic questions can 
be raised: (i) what are the limits of validity of such a non- 
equilibrium approach?, and (ii) are the mass fluctuation 
time scales observable? 
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FIG. 2: Asymptotic velocity distribution for different relax- 
ation time scales s of a Brownian particle, namely 0.1, 1 and 
10. The possible masses in the system are 1 and 10. The 
velocities are rescaled so that < v 2 >= 1. 
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IV. RELAXATION MECHANISMS 

In this section, we report numerical simulations of the 
random process Eq.QJ. The objective is twofold. First, 
we verify the theoretical prediction Eq.lJSJ in the adia- 
batic limit. Next, this allows to study systems which are 
beyond the range of validity of Eq.JSJl, namely systems 
where the separation of scales tn « t does not apply, 
thereby investigating the effects of competition between 
the relaxation to equilibrium and the fluctuating features 
of this equilibrium state. This program is achieved by 
considering three different relaxation characteristic time 
scales for the processes, namely s = ^f- = (0.1,1.0,10). 
Moreover, we first consider a paradigmatic case when the 
masses can switch between two different discrete values, 
each with equal probability: 

g(m) — —(6(m — mi) + S(m — 111,2)) (9) 



FIG. 3: Truncated Tsallis distribution gi(m) for a = 0.1, 
a — 1.0 and b — 00 in rescaled velocities, so that < |u| >= 1. 
The solid line is the Tsallis distribution -§y- — , , . The 

cut-off procedure gives more weight to the central range of 
the distribution (m g [1,3]), and decreases the importance 
of the peak and of the tail. The inset corresponds to these 
functions in log-normal scale, in the interval [0, 15]. 

If Eq.JSJ) applies, i.e. in the slow fluctuation limit, Eq.© 
leads to 



fs(v) = \\^-{^'ie-^ + V^e-^). (10) 
Z V Ztt 

Before focusing on the velocity distributions for the 
Brownian particles, let us stress that their average en- 
ergy depends on the speed of the fluctuation mechanism. 
Indeed, in the slow limit s << 1, the energy of the clus- 
ter converges very rapidly toward the equipartition value 
mi < v 2 >= e, where e is the average kinetic energy of 
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FIG. 4: Asymptotic solutions of the Langevin equation for 3 
time scales, where we fix the velocity scale < \v\ >= 1. The 
solid line is the theoretical distribution obtained from JSJ. 



FIG. 5: Stationary solution of the Langevin equation for 
s = 10 and truncated Tsallis distribution for a — 1. The 
dotted line is the Tsallis distribution — , n\ „,„ . 

7f (l + ( — ) 2 V 2 ) 2 



the bath and mj is the mass of the cluster at that time. 
This implies that the fluctuations of the measured veloc- 
ities are given by < v 2 >= 3(5^ + hi contrast, 
in the faster limit s > 1, the dependence of the charac- 
teristic time t fl ~ ^ can not be neglected. Indeed, this 
relation implies that particles with a smaller mass relax 
faster than particles with a larger mass. Consequently, 
the lighter particles should have a value close to their 
equipartition value rrii < v 2 >= e, while the heavier par- 
ticles should have an energy larger than their expected 
value. This property is verified (Fig. [JJ- 

In order to compare the velocity distributions for dif- 
ferent time scales r, and therefore at different energies, 
we rescale the velocities so that < v 2 > = 1. The results, 
as plotted in Fig. [21 confirm the theoretical predictions 
Ea. (|10fl and our description in the previous paragraph. 
Indeed, in the slow limit s = 0.1, the velocity distribution 
converges toward the leptokurtic distribution Ea. (|10|) . i.e. 
a distribution with a positive kurtosis and an overpop- 
ulated tail. In contrast, when the mass dependence of 
the fa has to be taken into account, most of the particles 
have velocities distributed like those of the light particles, 
namely are Maxwell-Boltzmann distributed. 

In the case of general Tsallis distributions, obtained 
through chi-squared distributed cluster masses, the prob- 
lem is more complex due to the continuum of masses in 
the system, and to the associated continuum of character- 
istic relaxation times. Moreover, the existence of extreme 
values for the masses may cause non-realistic numerical 
problems. Indeed, arbitrary small masses lead to arbi- 
trary high values of the velocities. In the following, we 
avoid this effect that is responsible for the power law tails 
observed in Tsallis distributions. This is justified by the 
fact that any physical system has a minimum size for 
its internal components. Similarly, we restrict the maxi- 
mum size of the clusters in order to avoid infinitely slow 
relaxation processes. These limitations are formalized by 



using the following truncation for the mass distributions, 
inspired by truncated Levy distributions p3j 



g(m; c) — k x 2 (to; c) if to > a and m < b 
g(m;c) = otherwise (11) 

where c is a parameter characterizing the chi-squared dis- 
tribution, k is a normalizing constant, and a < b are 
cut-off parameters. In the following, we use a = 0.01 
and b = 100.0 in the simulations of the Langevin equa- 
tion. As shown in Fig. this procedure is a natural 
way to smoothen the tail of the Tsallis velocity distri- 
bution while preserving its core. This method, that will 
be discussed further in a forthcoming paper, should be 
applicable to large variety of problems (like those men- 
tioned in the introduction) where extreme events have to 
be truncated for physical reasons. 

Numerical simulations of the Langevin equation for 
chi-squared mass distributions generalize in a straight- 
forward way the results obtained from the 2-level distri- 
bution Eq.@. Indeed, the faster the mass distribution 
fluctuates the larger deviations from the Tsallis distribu- 
tion are. Moreover, these deviations have a tendency to 
underpopulate the tail of the distribution, and therefore 
to avoid the realization of extreme values of the random 
process. 

It is also important to note the similarities between 
these asymptotic solutions of the Langevin equation and 
the truncated Tsallis distributions defined by Ea. (|llll . 
They highlight the flexibility of truncated Tsallis distri- 
butions in order to describe deviations from the Tsallis 
distributions, as shown in Fig. [S] 
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FIG. 6: Kurtosis of the distance distribution Eg. ll6L for 
two Bernouilli random walkers in one dimension with different 
walk lengths. 



V. DIFFUSION 

In this last section, we focus on the diffusive proper- 
ties of a Brownian particle with fluctuating mass. There- 
fore, we take into account the hydrodynamic time scale 
tni that is associated to the evolution of spatial inhomo- 
geneities. It is well known 24], when the mass of the 
Brownian particle is constant in time, that separation of 
time scales tn » tn is required in order to derive the 
diffusion equation: 



dtn(x;t) — Dd^nfat) 



(12) 



where D is fixed by the Einstein formula. 

When the mass of the Brownian particle fluctuates, 
however, we have shown above that there is an additional 
time scale r in the dynamics. As first approximation, we 
restrict the scope to the limit of very slow fluctuations 
(r >> tn), that is more restrictive than the limit dis- 
cussed in the previous section (t >> In that case, 
it is possible to show [2j| that the Chapman-Enskog pro- 
cedure leads to Eq.{7|), where D is now time-dependent, 
i.e. a random variable that is a function of the mass of 
the Brownian particle. 

In the following, we investigate the process associated 
to Eq.lQ with fluctuating diffusion coefficient. To do so, 
we simplify the analysis by considering a one-dimensional 
discrete time random walk, where the jump probabilities 
may fluctuate in time j26(. Namely, the walker located 
at x performs at each time step a jump of length I, with 
Bernouilli probabilities: 



P(k)\ l = ±[6(k,l) 



■S(k,-l)} 



(13) 



The quantity / fluctuates between two integer values 
I a < Ib that correspond to an heavy/light state for the 
Brownian particle. When I a — Ib, it is easy to show that 
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FIG. 7: Kurtosis k of the distance distribution, as a function 
of time, for the Bernouilli random walker with fluctuating 
mass (see main text). In (a), we focus on the limiting cases 
p = 0.0 and p — 1.0 that converge toward the theoretical 
prediction, Ea. l|16|l . « = 1.08 and kg = (Gaussian) re- 
spectively. These asymptotic values are represented by solid 
lines. In (b), we focus on the intermediate cases p = 0.01 
and p = 0.0001. The state associated to Eg. 11611 is stationary 
during a long time extent, and is followed by a convergence 
toward the Gaussian. 



the first distance moments rrii 
behave like: 



-< x % > asymptotically 



m% = Dt 
m 4 = 3D 2 t 2 



(14) 



where the diffusion coefficient D = l 2 A 



Relation 113)1 

implies that the kurtosis of the distance distribution 
K _ mi _ g vanishes asymptotically, as required by the 

m 2 

central-limit theorem. 

Another simple limit consists in a system with I a ^ Ib, 
and where the mass fluctuation process is infinitely slow. 
Consequently, the particles do not change mass and the 
system is composed of two species that diffuse differently. 
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In that case, the first distance moments read: 



rri2 
777,4 



\{Da 



D B )t 
D%)t 2 



(15) 



where the diffusion coefficients are Da = l\ and Db = 
l 2 B . This relation implies that the asymptotic kurtosis is 
equal to: 



k = 6 



D%) 



{Da + DbY* 



(16) 



In Fig© we plot n{r) = - 3, where r = The 

figure shows that k > 0, i.e. the distribution is char- 
acterised by an overpopulated tail, except in the usual 
limit r = 1. This result is expected, as a system with 
r 7^ 1 is composed of two species that explore the space 
at different speeds. Let us stress that despite this anoma- 
lous position distribution, the diffusion is standard, i.e. 
< x 2 >~ t (see Eq[T3J). 

In order to study intermediate situations, we have per- 
formed computer simulations of the random walk. The 
system is composed of 50000 walkers, that are initially 
located at x — and randomly divided in the species 
A/B. The mass fluctuations are uncorrelated and occur 
with probability p G [0, 1] at each time step, i.e. the 
characteristic time of the fluctuations is r ~ pT 1 . In the 
simulations, we have used I a = 1 and Ib — 2, so that 
the prediction for the kurtosis in the limit r — > 00 is 
kab = 1.08. This prediction (FigEJi) is verified by simu- 
lations in a system where the mass of a walker is constant 
in time. In the limit r — 1, i.e. the walker changes mass 
at each time step, the kurtosis converges rapidly toward 
the gaussian value kq — 0. In contrast, for higher finite 
values of r (FigUJa), one observes a crossover between the 
two asymptotic behaviours kab and kq. The stability of 
the state kab is longer and longer for increasing values of 
t. This quasi-stationarity j2]j originates from the follow- 
ing reason. Over a long time period T, with T » r, the 
particles have suffered so many changes from state A to B 
and forth, that their asymptotic dynamics is equivalent 
to that of a random walk with 2 possible jump lengths 
I A and Ib, with probabilities: 

P(k) = ^[S(k,l A )+S(k,~lA)+S(k,l B )+S(k,~lB)} (17) 

Consequently, its asymptotic dynamics is subject to the 
classical central-limit theorem, and the position distribu- 
tion is a Gaussian for T >> t. 



the simplicity is related to underlying considerations on 
quite various physical (or other) systems in which some 
"mass" is evolving with time, sometimes stochastically, 
as mentioned in the introduction of this paper. 

For the sake of generality, we have studied systems 
with arbitrary mass distribution and time scale of the 
mass fluctuation. Let us stress that the mass statistics 
has been chosen a priori, and that a more detailed study 
requires dynamical treatment of the fluctuating mass by 
Langevin equation. This additional modelling depends 
on the nature of the considered problem. The velocity 
distribution of the Brownian particle has been studied by 
performing simulations that highlight the important role 
of mass fluctuation time scale. In the case of chi-squared 
mass distribution, it is shown that truncated Tsallis dis- 
tributions seem to describe in a relevant way deviations 
from the Tsallis statistics. Such distributions should be 
applicable to a large variety of problems where extreme 
events have to be truncated for physical reasons, e.g. fi- 
nite size effects, - when there is no infinity! Among these 
possible applications, let us note their occurrence in air- 
line disasters statistics 28] . 

We have also studied the diffusive properties of the 
Brownian particle. In the limit of slow mass fluctuation 
times, the particle motion is modelled by a random walk 
with time-dependent jump probabilities. Moreover, for 
the sake of clarity, we restrict the scope to a dichotomous 
mass distribution. This modelling is a simplification of 
the complete problem, that has to be justified by a com- 
plete analysis starting from the Fokker-Planck equation 
itself [28|. Nonetheless, despite its apparent simplicity, 
the random walk analysis shows non-trivial behaviours, 
namely the system is characterised by standard diffusion, 
associated with non-Gaussian scaling distributions and 
quasi-stationary solutions. These features originate from 
very general mechanisms that suggest their relevance in 
various systems with fluctuating "mass" parameters. 

In conclusion, we have generalized Brownian motion 
to the case of fluctuating mass systems, and distinguish 
them from systems in which there are local temperature 
fluctuations. We have found that the velocity of such 
a particle may exhibit various anomalous distributions, 
including Tsallis distributions and truncated Tsallis dis- 
tributions. The study of a one dimensional random walk 
also indicates anomalies in the kurtosis of the position 
distributions which might indicate physical biases in var- 
ious processes as those recalled in the introduction. This 
suggests to look for the value of higher order distribution 
moments and on their time evolution for understanding 
properties of the non-equillibrium systems. 
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